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Abstract 



We present a new method to evaluate with high precision leading isospin breaking effects due to the 
small mass diflFerence between the up and down quarks using lattice QCD. Our proposal is applicable 
in principle to any hadronic observable which can be computed on the lattice. It is based on the 
expansion of the path-integral in powers of the small parameter rrid — ruu. In this paper, we apply 
this method to compute the leading isospin breaking effects for several physical quantities of interest: 
the kaon meson masses, the kaon decay constant, the form factors of semileptonic Kis decays and the 
neutron-proton mass splitting. 



1 Introduction 



Isospin symmetry SU{2)y^ is an almost exact property of strong interactions as described by the QCD 
Lagrangian. This happens because the difference between the up and down quark masses is much smaher 
than the QCD scale, {rrid — mu)/^QCD <C 1, and it remains true also when electromagnetic interactions 
are switched on, because isospin breaking effects due to the different quark electric charges {cu 7^ e^) are 
suppressed by the electromagnetic coupling constant, Ofem ^ 1/137. For these reasons most of theoretical 
predictions of several physical quantities assume isospin symmetry, i.e. the masses of the up and down 
quarks are taken equal and electromagnetic effects are neglected. 

Nowadays, with the increasing precision of the experimental determinations of many physical quan- 
tities, and in some cases with the improvement of the theoretical predictions, the control over isospin 
breaking effects is becoming phenomenologically relevant. This is the case, for example, of the form fac- 
tors parametrizing Kis decays. Isospin breaking effects are important also for hadron spectroscopy, for 
the meson decay constants, for the tt-tt scattering length, for the quark condensate and for many other 
quantities. 

In the past, isospin breaking effects due to the light quarks mass difference (in the following referred to 
as QIB effects for QCD isospin breaking effects) have been accommodated within the chiral perturbation 
theory (xpt) framework (see 1-6 for a largely incomplete list of references on the subject), while several 
attempts to compute electromagnetic effects for the hadron spectroscopy in lattice QCD have been pre- 
sented [7[|1Q]. It is very difficult to take into account in numerical simulations QCD isospin breaking (see 



refs. 11-15 for a selection of previous lattice works on the subject) because the effect is in general rather 
small and comparable with the errors in the determination of, say, the meson masses or decay constants. 
Furthermore, in order to perform unitary dynamical simulations of two light quarks of different mass the 
single quark determinant must be positive and this happens only in the case of lattice discretizations of 
the fermion action that are very expensive from the numerical point of view. 

In this paper we present a new method to compute the leading QIB effects with high precision. The 
method is based on the expansion of the lattice path-integral in powers of the small parameter rrid — ruu 
and is applicable in principle to any hadronic observable which can be computed on the lattice. As a first 
application, and to show that the method works, we have applied it to the computation of the leading QIB 
effects for several physical quantities of interest: the kaon meson masses, the kaon decay constants, the 
form factors of semileptonic K^s decays and the neutron-proton mass splitting. In the future we plan to 
apply the method to other physical quantities, to include QED corrections and to try also the calculation 
of next-to-leading corrections such as the 7r+-7r^ mass difference. 

The main physical results of this work are 

[m,-mJ^^^(M5,2GeF) = 2.35(8)(24) MeV x 



Fk/F^ 



6.05 X 103 MeV^ 

= -0.0039(3)(2) X ^f^^ , , 

^ ^ 6.05 X 103 MeV^ 

.2 iQCD 



[M^-M^]^^^ = 2.8(6)(3)MeV x M^_J^±J , (1) 
^ ^ ^ 6.05 X 103 MeV^ ^ ^ 

and the QIB corrections to semileptonic decay rates discussed in sec. jT] In previous equations the lattice 
error (the one in the first parenthesis) has been obtained with a rather modest statistics, ~ 150 gauge field 
configurations, and can be reduced in the future. Most of the systematic error comes from the ambiguity 
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in the definition of the electromagnetic corrections (see uncertainties in the second parenthesis). In this 
study we have extracted the rrid — ruu mass difference in QCD by using as an external input the QIB effect 
in the kaon masses, [M^o - M]^+]^^^. This quantity requires in turn a definition of the electromagnetic 
isospin breaking corrections that we have taken from other calculations, as explained below in the paper, 
and that will be treated more properly in the future. The result for the neutron-proton mass difference 
has been obtained at fixed lattice spacing and will be improved in a separate publication. 

Next-to-leading QIB corrections, i.e. second or higher orders in the (m^ — rrid) expansion, have not 
been calculated in the present paper. The associated effects are estimated to be negligible at the current 
level of (both theoretical and experimental) precision on flavour physics observables. Indeed, on the basis 
of dimensional arguments, higher orders corrections are expected to be suppressed by additional powers 
of the small expansion parameter {rud — Tnu)/^QCD- The actual numerical size of the leading order 
contributions presented in eq. ([T]) is consistent with this power counting expectation. Clearly, we cannot 
exclude the existence of specific observable for which higher order QIB effects may be larger than expected. 
Nevertheless if necessary, the method suggested in the present paper can be extended to evaluate higher 
order contributions as discussed, for example, in sec. [3] for the 7r+-7r^ mass splitting. 

In the case of and F^, we quote also the first order derivatives with respect to rrid — m^. 



[md-mu]^c^{MS,2GeV) 



K+1 _ or7/QAx/in3 



2.57(8) X lO'^ MeV, 



--IL^^'' = 3.3(3) X 10-^ MeV-^ (2) 

so that one can use his/her preferred value of the up-down mass difference to get the physical QIB effect. 

The paper is organized as follows: in sec. [2] we introduce the method and set our notation. In sec.|3]we 
derive the correlation functions needed to extract isospin breaking effects whose calculation is discussed 
in detail in the remaining sections. In particular, in sec. [4] we study kaon two point correlation functions 
and in sec. [5] we discuss electromagnetic isospin breaking effects and extract rrid — and F^o — Fx+' In 
sec.|6]we calculate the neutron-proton mass difference while in sec.[7|we discuss the calculation of isospin 
breaking effects for the K^s form factors. Our conclusions and outlooks are given in sec. [Sj Some technical 
details are discussed in the appendices. 



2 Description of the method 

In this section we present the basic ingredients of our method, which is simply based on a perturbative 
expansion in the small parameter {md — mu)/^QCD- Let us start by considering the evaluation of a generic 
euclidean correlation function (O) used to extract information about physical quantities as masses, decay 
constants, form factors etc., 

where Dcj) represents synthetically the full functional integration measure of the theory. By neglecting for 
the moment electromagnetic corrections and possible isospin breaking terms that may arise because of 
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lattice artifacts with particular choices of the lattice fermion action, we can write the Lagrangian density 
as a term which is SU{2)y symmetric plus a term which violates the isospin symmetry. 

- ruu + Tad / - T 7X Tad — / _ t 7X 
= Ckin H ^ [uu + dd) [uu - dd) 

= Ckin + raud qq - ^raud qr^q 

= Cq- Amud t , (4) 

where = (u^d), niud = + ^u)/2 and ArUud = (^d ~ By expanding at first order the 

exponential of the action, S = with respect to ArUud we obtain 

^ JD4>0{1 + Am^d S) e-^° ^ (0)o + Am„rf {OS)o 
JD^ {1 + Amw S) e-So l + Am^ (^)o 

= {0)o + Amud{OS)o, (5) 

where (•)o represent the vacuum expectation value in the isospin symmetric theory and S is the isospin 
breaking term. 



The correction in the denominator vanishes, (S)o = 0, because of isospin symmetry. Concerning the Wick 
contractions of the correlation functions {OS)o, isospin symmetry makes also to vanish some fermionic 
disconnected contributions of the form 

Aniud { [fermionic Wick contractions of O] x tr[6'] )o = , (7) 

since these are proportional to the trace of the flavour matrix t^. We can now describe a general recipe to 
be used in order to compute leading QIB effects on the lattice: 

• consider a given correlation function in the full theory, i.e. with ^ md^ and for each gauge 
configuration draw all the fermionic Wick contractions; 

• expand the up and down quark propagators with respect to Aniud according to 

Gu{xi,X2) = G(,{xi,X2) + Amud ^Gi{xi,y) Gi{y,X2) ^ , 

y 

Gd{xi,X2) = Gg{xi,X2) - Amud ^ Gg{xi,y) Gg{y, ^2) H ; (8) 



3 



• retain the terms linear in ArUud and compute the corresponding diagrams (or fermionic Wick con- 
tractions). 

In the fohowing sections we shah discuss in detail how to extract physical information from the resulting 
correlation functions. To this end we need to set the notation we are going to use in drawing diagrams. 
Eqs. P can be represented diagrammatically as follows 



d 




where here and in the following the up quark line in the full theory is drawn in light blue color while the 
down quark line in green. All the black lines refer to G^, the propagator with the symmetric mass rriud 
in the isospin symmetric theory whose square, entering eqs. ([8|, can be easily calculated on the lattice by 
using Gi itself as the source vector of a new inversion. 

The insertion of the scalar density is represented by a cross according to 
y > , = Ge{x - y) = {i{x)i{y)) , 

Z Z 

with £ either u or d. Here and in the following the superscript L stays for bare lattice quantity. In particular 
we have 

Arriud = Zatti ^rn^d (H) 

where Z/\rn is scale and scheme dependent while the combination Arriud ^(^)^(^) renormalization 
group invariant. According to eq. ([s]), a generic correlator {OS)o with a single insertion of the isospin 
breaking term can be obtained from (O), a correlator in the full theory, as the derivative of the latter 
with respect to Arriud evaluated at Arriud = 0. It follows that, by working within a mass independent 
renormalization scheme, {OS)o is finite provided that O and S have been separately renormalized. By 
iterating the previous argument it can be easily understood that the connected parts of correlators with 
multiple insertions of the renormalized operator S are finite. 

In this paper we have applied the method discussed above by using the so called Twisted Mass lattice 
discretization of the QCD action. This choice has advantages and drawbacks. The big advantage is auto- 
matic 0(a) improvement. The drawback is the breaking of isospin symmetry at finite lattice spacing even 
with Arriud = 0. The associated O(a^) cutoff effects are eliminated by performing continuum extrapola- 
tions. Since our method is general and can be applied with any lattice regularization of the quark action 
(e.g. Wilson, Overlap, etc.), in the main body of the paper we illustrate our results without entering into 
the specific details of the Twisted Mass fermion action that we discuss in Appendix [Aj 

We close this section by explaining the notation used in the following to express and calculate variations 
of correlation functions and matrix elements. To this end it is useful to introduce the following operators 
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acting as absolute variations 



A/O = 0(f ) - 0(i} , 
AbO = 0{{) - 0{h) , 

^^^A,0 + A,0^0if)-0ib)^ ^^^^ 



and the corresponding relative variations 

AfO A,0 AO 

Here and 0{b) are quantities calculated with the light quark propagators at first order in Arriudi 

while 0{i) is the corresponding quantity calculated in the unperturbed isospin symmetric theory. The 
labels / (forward), b (backward) and i (iso-symmetric) are generic. More precisely, / will stand for 

f = {d, K^, n, K^t:- ] (14) 

when we shall discuss in turn quark masses, kaon masses and decay constants, the neutron and proton 
masses, the semileptonic decay D Klu and the semileptonic decay K irlu. Correspondingly b and i 
will represent 

b = {u, K^, p, D^K"', K^TT^ } , 

i = {£, K, N, DK, Ktt } , (15) 

where TV stands for nucleon. 



3 Correlation functions at first order 

In this section we shall derive the correlation functions that need to be calculated in order to extract 
the leading QIB corrections to meson masses and decay constants, nucleon masses, and to the form 
factors parametrizing semileptonic meson decays. In particular, we shall consider the following two point 
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correlation functions 



X 

X 
X 

CKOKo{t,p) = J2 e-'^-^{dj,s{x) S75d(0) ) , 



X 



(X) 
(X) 



edef- 



1 ±7'- 



Udiu^C^^df) 



1 ± 7° 

^def 7, dd (d^ C-f5 Uf ) 



(0)), 

(0)), 



and the following three point correlation functions 

x,y 

C^^Koit-^PD,PK) = ^e-'P--^e-'P-<^-y^d^''c{y,T/2) crs{x,t) s^^'diQ)) , 

x,y 

C^o^-{t;PK,P.) = ^e-'P--^e-'P-<^-y\d^''s{y,T/2) s^>^u{x,t) u^HiQ)) , 

x,y 



(16) 



(17) 



In previous expressions C denotes the charge conjugation matrix while €abc the totally antisymmetric 
tensor in color space. 

A first trivial observation comes from eq. ^ telling us that all the quantities that do not involve a 
light valence quark propagator do not get corrected at first order in Amud- This is the case for example of 
heavy-heavy and heavy-strange meson masses and decay constants, etc. Pion masses and decay constants 
too do not get corrected at first order. This can be shown diagrammatically for the charged pions two 
point function 



0^0 



(18) 



and for the connected diagrams entering neutral pion two point function C-j^o-j^o{t)^ 



2 <^^^ + 0{Am„d? , 



2 <^^^ + OiAniudf , 



(19) 



It is easy to show that the first order corrections cancel also for the disconnected diagrams contributing to 
C^o^o {t) in the full theory, a known result that can be understood in terms of isospin quantum numbers. 
The Wigner-Eckart reduced matrix element of the operator S between pion states is indeed zero for G- 
parity, 



(^||5|K) = (l,/3||l,0||l,/3)=0. 



(20) 



This is certainly not the case at second order where the relevant 0{Amud)'^ diagrams are 



(21) 



For flavoured mesons first order corrections to masses and decay constants are instead different from 
zero. Here we discuss the case of strange particles but the discussion proceeds unchanged if the strange is 
replaced with a charm or a bottom quark. The QIB correction to the two point correlation functions of 
the strange mesons are 



(22) 



In the diagrams above and in the following the strange quark line is red. Note that the correction to the 
neutral mesons is equal in magnitude to that to the charged particles, AfCKK{t) = AbCxKit) = ACxKit)- 



We now consider first order corrections in the case of nucleon masses. The neutron-proton mass 
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difference can be extracted at first order in Arriud by the diagrammatic analysis of C^^{t) and Cpp{t)^ 



1=1 F=I-4F=I F=l f=) 



(><) - [><3 = 4|><) - - 1><3 



(23) 



As usual this is obtained by expanding all the light quark propagators appearing into the correlation 
functions. Also in this case we find AfCNN{t) = AijCNN{t) = AC n Nit)- 

Concerning the form factors parametrizing semileptonic decays, we start by considering a charmed 
meson decaying into a strange meson. The discussion would proceed along the same lines in the cases of 
B ^ D transitions. The charm quark line is drawn in yellow. We get 




+ 0(Am„d)' , 



0(Am„d)2 . 



(24) 



As for the correlation functions analyzed above the correction is equal in magnitude between the two 
processes, because the weak flavour changing current does not contain a light quark field. From the previous 



two equations one can extract /+°^^(^^) — /+^^ (^^)- At present the experimental and theoretical 
uncertainties on this quantity are such that QIB effects can be safely neglected but we have chosen to 
discuss these processes first because of their simplicity. Indeed, the phenomenologically relevant K ^ tt 
case involves disconnected diagram^ that complicate the analysis. 

The expansion of C^o^_ (t) is given by 




0(Am„d)2 . 



(25) 



■"^This is not in contradiction with eq. ^ since disconnected diagrams arise by making the fermionic Wick contractions of 
the observable O and not from TrlS*]. 
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The correction in the tt^I^v case is obtained from the correlation function C^^^Q(t) in the fuh 

theory, whose disconnected diagrams survive at first order in ArUud^ 



s>0 



O 



O 




(26) 



The correction to 7t^I~^u is not just equal to the tt one. From eqs. (26) and (25), one 



can extract f^^^ {q^) — ^°(^^)- At = this quantity has been estimated by using the measured 
form factors and chiral perturbation theory [16^ , 



"(0) 



-0.029 ±0.004, 



(27) 



and found of the same size of the deviation from unity of the form factor at q'^ = 0, being ff^{0) 
/(f ^(0) = 1 the value at the SU{3)v symmetric point {Mk = M^). 



4 Kaon masses and decay constants 



In this section we discuss in detail the strategy used to derive the isospin corrections to the kaon masses 
and decay constants. To this end we start by considering the Euclidean correlation functions of eqs. ( [22| ) 
both in the full theory and in the isospin symmetric one. The spectral decomposition of Cxo^o (the 
analysis of Ck+k- proceeds along similar lines) is 



c. 



e-'P-^d^,s{x, t) S75rf(0)) = 



{0|d75g(0)|n^) (n^ 1^75^(0) |0) 
2E^ 



(28) 



where the dots represent sub leading exponentials and where |n^) and are the states and the eigenvalues 
of the perturbed theory corresponding respectively to |n) and En in the isospin symmetric unperturbed 
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theory. These are related at first order in perturbation theory with respect to ArUud according to 



= = \K) + \AK) + OiAm^df , 



Eko =E^ = Ek + AEk + 0(Am„d)2 . 



(29) 



Explicit expressions for AEk and \AK) are derived in Appendix [B| By substituting eqs. (29) into eq. (28) 
and by recalling the diagrammatic analysis of eqs. (22) we obtain 



^-EKt 



2Ek 

^-EKt 



2Ek 



A{GI/2Ek) 
Gl/2Ek 



tAE, 



K 



(30) 



where AGk = (0| 575^(0) |AK). Note that the insertion of the QIB term jC constitutes a flavour diagonal 
perturbation and that, consequently, the kaons are the lightest states contributing both to Ckk{p^ t) and to 
ACkk{p-, t). The analysis would be considerably more complicated in the case of a perturbation (typically 
insertions of the weak hamiltonian) opening a decay channel for the kaons because the physical information 
would be hidden into sub-leading exponential terms. 



In our case, by studying the ratio of the two correlators of eqs. (30) 



(31) 



it is possible to extract the leading QIB corrections to kaon energies and decay constants. Indeed AEk 
appears directly in the previous equation as the "slope" with respect to t whereas SFk can be extracted 
from the "intercept" according to 



Fk {rus ^mud)j^ 



K 



SFi 



Am. 



K 



ud 



rris + rriud 



SGk - 2SMk . 



(32) 



On a lattice of finite time extent T with quark fields satisfying anti-periodic boundary conditions along 
the time direction and given our choice of the kaon source and sink operators, the pseudoscalar densities. 



eq. (31) has to be modified according to 



8CKK{m = ^' 



2E 



K 



• AEK{t - T/2) tanh [EK{t - T/2)] 



(33) 



10 



0.6 [ — ^ ' ' ' ' ^ ' ' ' ' ^ ' ' ' ' ^ ' ' ' ' 1 25 [ — r 



r 




Figure 1: Left panel: extractio n of meson energies from the effective mass of CxKiv^t)- Right panel: fits of 
^CKKip^t) / ci^^ud according to eq. jssj: as it can be seen numerical data follow theoretical expectations. The data 
correspond to (3 = 3.9, am^^ = 0.0064, am J = 0.0177 (see Appendix [a). 




Figure 2: Top-left panel: the effective energies extracted from CxKip^t) are compared with the expected continuum-like 
dispersion relation E'^(p) = p"^ -\- (no fit). Top-right panel: the corrections to the effective energies are compared with 
the theoretical expectation AEk(p) = MkAMk/Ek(p) (no fit). Bottom-left panel: Fk extracted from CxKip^t) (dark red 
points) and from C^j^(p,t) (dark blue points, see eq. ( |34| ) for different values of p; the solid line correspond to the value 

extracted at p = 0. Bottom-right panel: SFk /aAm^^ extracted from SCKK(p,t) (dark red points) and from 6C^^(p,t) 
(dark blue points) for different values of p; the solid line correspond to the value extracted at p = 0. The data correspond to 
(3 = 3.9, am^^ = 0.0064, amf = 0.0177 (see Appendix[A|). 



As can be seen from Figure [T] SCkk{p^ t) is determined with high precision, given the strong statistical 
correlation existing between the numerator and the denominator of the ratio in eq. (31). A consistency 
check of our procedure consists in verifying the dispersion relation Ej^{p) = + M|^ and in comparing the 
variation AEk{p) against its expected behaviour AEk{p) = Mk AMk /_Ek{p)' Excellent agreement is 
found between numerical data and the theoretical curves shown in Figured both for E'^{p), top- left panel, 
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and AE'^{p)^ top-right panel. In the bottom panels of Figure [2] we also show that two different definitions 
of 5Fk (blue and red points) extracted from correlators at several p- values give consistent results. The 
second definition of Fk and of 5Fk has been obtained by considering the correlation function between the 
pseudoscalar density and the axial vector current 

CfAt) = Y.^h'l'^<^)-si,m) = + (34) 



and its correction at first order in Arriud- 

We shall now discuss the chiral and continuum extrapolations. Concerning chiral extrapolation, it is 
useful to consider the correction to the meson mass square because this is a finite quantity in the chiral 
limit. The chiral formulae for AM^ and SFk have been obtained long ago, in ref. within the unitary 
SU{3)l X SU{3)r effective theory (iV/ = 2 + 1), 

AMI _ , 2^^^^^ , 



= 5o| 



+(ms + ■mud)—^{2Ll - Lg) + (m^ + 2m„rf)— 2-(2Le - L^j 
^0 ^0 



4)} , 

SFk _ „ f 4L| _ _l hk_ _ l^r, - 



where Bq, Fq and are low energy constants while 

M2 = 2BoTO„d , 

M|- = Boiniud + nis) , 
M2 =2Bo(m„d + 2m,)/3, 



(35) 



/U-2 

/ip = ^^;^ln(M|,/M'), P = {7r,v,K}. (36) 



In view of the poor convergence properties of SU{3)l x SU{3)r effective theory and also to cope with 
the fact that our results have been obtained by quenching the strange quark {Nf = 2), we have chosen 
to fit our data with the formulae resulting by expanding the r.h.s. of eqs. ( [35| ) with respect to rriud/'^s 
(see ref. [TT^ for a detailed and quantitative discussion of this point). This procedure is justified when the 
average light quark mass is sufficiently small compared to both Aqcd and the (valence) strange quark 
mass, as it appears to be the case by looking at our data shown in Figure [3j 

More precisely, in the left panel of Figure [3] we show the combined chiral and continuum extrapolation 
of AM]^/ Arriud' Our results, obtained at four different lattice spacings for several values of the average 
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Figure 3: Left panel: combined chiral and continuum extrapolations of AM^ / Arriud- Right panel: combined chiral and 
continuum extrapolations of SFx / Arriud- Black points correspond to the coarser lattice spacing, a = 0.098 fm, dark magenta 
points correspond to a = 0.085 fm, green points to a = 0.067 fm and blue points to a = 0.054 fm. Red lines are the results 
of the continuum extrapolations. 



light quark mass (see Appendix [A|, do not show a visible dependence with respect to ruud that can 
be quantified within the quoted errors. These errors combine in quadrature the statistical error and the 
systematic one coming from the uncertainties on the lattice spacing and from the renormalization constants 
(see Appendix [A|. We have consequently fit the numerical data to the following functional form 



Am 



ud . 



Am. 



ud 



QCD 



(37) 



where [AM'^/Arriud]^^^ is a constant thus representing the physical value of this quantity in the contin- 
uum limit and at the physical light quark mass. We have also performed a fit of the numerical data by 
adding to the previous expression a linear and a (chiral) logarithmic term. In this case we have checked 
that the fitted coefficient of the chiral logarithm is compatible within its (large) error with the value 
—Bq/{S7t'^Fq) resulting from the expansion of eqs. (35) with respect to rriud/'^s' The two fits give com- 
patible results for [AM]^ / Amud\^^^ within the associated errors and we consider the difference between 
the central values as an estimate of the chiral extrapolation systematics, an error of about 3% that we 
have added in quadrature to the lattice uncertainty. 

In the case of SFk/ Arriud^ right panel of Figure [sj the dependence upon rriud is significant within 
quoted errors (again combining in quadrature statistical and systematical ones) and we have included in 
the fitting function the leading and next-to- leading terms appearing in eqs. (35) expanded in powers of 
rriud/'^s plus a lattice artifact term, i.e. 



SF 



K 



Am,, 



{rriud^a) 



SF 



K 



Am,, 



QCD 



CfO^ 



rriud 



m 



QCD\ 
id J 



B2 rriud log 



rriud 

QCD 
ud , 



m 



(38) 



The systematics associated to this extrapolation has been estimated by replacing the (chiral) logarithmic 
term appearing into the previous expression with a quadratic term and it has been found of the order of 
5%. 



The value m! 



QCD 

ud 



m'^^''{MS,2GeV) 



3.6(2) MeV has been taken from refs. 18 19 . The fitted 
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values 



QCD 

(MS, 2GeV) = 2.57(8) x 10^ MeV , 

QCD 

(MS, 2GeV) = 3.3(3) x 10"^ MeV"^ , (39) 



will be used in the next section to obtain our physical results for Am^^^ = Am^^^ {MS,2GeV) and 



5 Electromagnetic corrections and ma — niu 

In this section, by using as input the experimental determination of M'^o — , we shall determine 

When comparing the theoretical predictions with the experimental numbers we cannot neglect the 
isospin breaking corrections induced by electromagnetic interactions. In the literature, it has become 
popular to separate electromagnetic and QIB effects and to give separately the theoretical value of these 
two contributions. Different calculations of the electromagnetic corrections, performed on the lattice or 
within other non-perturbative approaches, are also often compared |18 . However the separation of the 
electromagnetic and strong QIB effects is ambiguous [2Q|[2T] , i.e. the different contributions depend on the 
definition by which they are separated, whereas they do not correspond to any physical observable because 
of ultraviolet divergences. 

If we work at first order in the QED coupling constant and Arriud and neglect terms of 0{aemAmud)^ 
the relevant diagrams entering the difference of kaon two point functions are 



Aruud 



Arriud 




The electromagnetic corrections to CKK{t) are logarithmically divergent, corresponding to the renormal- 
ization of the quark masse^ and the divergent part has two components: one proportional to qq and 
the other proportional to qr^q. Alternatively one can treat electromagnetism on the lattice to all orders 
by exploiting the QED non compact formulation 7 but, in any case, two independent renormalization 
conditions have to be imposed in order to fix the counter-terms and separately renormalize the up and 
down quark masses. This can be achieved by matching the physical masses of the charged and neutral 
kaons of the present example (the mass of the strange quark could be eventually fixed by the mass of the 
baryon). Having extracted the light quark masses, i.e. rriud and Arriud^ one can predict F^o —Fx+^ the 
proton and neutron masses and all the other observables by including "physical" isospin breaking effects. 

In this paper we only consider the QCD corrections because we want to show that our method works for 
this part (we shall present a proposal for the 0{aem) corrections in a separate paper). This is equivalent to 

^For simplicity we neglect in the discussion the renormalization of the meson sources/sinks, which in any case are finite 
if we use vector or axial currents. 



14 



say that we follow the common procedure of separating the two isospin breaking contributions by switching 
off electromagnetism. Obviously the attempt to use physical quantities to fix Am^^^ fails, since there 
are no data with electromagnetic interactions switched off and, for this reason, we shall use the definition 
and determination of the electromagnetic corrections worked out by other groups. 



According to Dashen's theorem [22], electromagnetic corrections are the same in the chiral limit for 
f|^+ and M^o — while, as discussed in the previous sections, pion masses are not affected by 
first order QCD corrections. Beyond the chiral limit it is customary |T8| to parametrize violations to the 
Dashen's theorem in terms of small parameters and, concerning M'^o 



Ml, 



M|.+ , we have 



[M|o - Ml,] 



2 ^QCD 



[M|o-M^,]^--(l + e,) [Mlo-Ml,r\ 



(41) 



where we have neglected QCD contributions of the second order 0{Am'^^) in the pion mass difference. 
By using chiral perturbation theory and results from lattice QCD calculations of the electromagnetic 
corrections [7 10 , ref. 18 estimates 



0.7(5), 



[Mlo - M|+]^^^ = 6.05(63) X 10^ MeV^ . 



(42) 



On the one hand, it is important to realize that to fix a value for the electromagnetic term is equivalent 
to define a prescription to separate (first order) QCD and QED isospin breaking effects and that for this 
reason we might ignore the error on this quantity. On the other hand, one may wonder whether the 
different results combined in ref. flS^ have really been obtained by using the same prescription. For this 
reason, in the following, we shall quote our results by considering the error on as a way of taking largely 
into account this "scheme uncertainty" (see ref. [20] for a detailed discussion of this point). 



By using eqs. (42) and the results for [AM'l/ Arriud]^^^ and [SFx / Arriud]^^^ given in eqs. (39), we 
get the following results 



[rrid 



iQCD 



{MS, 2GeV) 



2.35(8) (24) MeV 



[M 



\QCD 



6.05 X 103 MeV^ 



QCD 



0.0078(7)(4) 



]QCD 



K+1 



6.05 X 103 MeV^ 



(43) 



where the first error comes from our calculation and combines in quadrature statistics and systematics 
while the second comes from the uncertainty on s^. 

At first order in ArUudi thanks to the fact that pions don't get corrections and that and get 
opposite corrections, we have 



QCD 



= -0.0039(3)(2) 



[M^ 



Ml, 



\QCD 



6.05 X 103 MeV^ 



(44) 
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a value that is higher than the estimate obtained in ref. 6 by using chiral perturbation theory, namely 



L FkIf^ 



1 



-0.0022(6) . 



(45) 



6 Nucleon masses 

Having determined Am^^^ , we can now predict the QCD isospin breaking corrections on other observ- 
ables, as already done in the previous section for the kaon decay constants. In this section we calculate 
the difference between the masses of the neutron and of the proton. We consider the correlation functions 
C^^(t) and Cpp{t) at zero momentum (see eqs. (16)) and, in order to decrease statistical errors, we extract 
nucleon masses from the combinations 

Cppit) = C+it)-C-^iT-t). (46) 
The quark fields entering the sink interpolating operators have been "Gaussian smeared" according to 

^smeared ^ (I agH^H , 
3 

H{x,y) = ^^Ui{x,t)5^^^-i^Uj{x-i,t)S^^^^i , (47) 



with the parameters ag and Ug fixed at the values optimized in ref. 23 where the same gauge configurations 
of this study have been used. 

The extraction of physical informations from nucleon euclidean two point functions proceeds along the 



same lines described in detail in the case of the kaons. By using the diagrammatic analysis of eqs. (23) we 
have 



Cat AT (t) 



(48) 



and 



SCNN{t) 



S ^ B ^j^Sira 

3 



j^^- 













SWn - tAMN 



(49) 
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Figure 4: Left panel: effective mass plots of nucleon correlation functions CiVA^(t). Right panel: Correlation functions 
5C N N {t) / a/S.m^^. The data are at fixed lattice spacing a = 0.085 fm for different values of rriud (see Appendix [a|> . 




20 30 40 

m,d (MeV) 



Figure 5: Chiral extrapolation of AMjy / Arriud at fixed lattice spacing a = 0.085 fm. 



where the dots represent sub-leading exponentials contributing to the correlation functions. By extracting 
the slope in t of SCNN{t)^ we can determine AMn = {Mn — Mp)/2. In the left panel of Figure [4] we show 
the effective masses M^^ (t) as extracted from the correlation functions CNN{t) at fixed lattice spacing 
a = 0.085 fm for different values of mud- In the right panel of the same figure we show the corresponding 
correlation functions 6CNN{t) /aAm^^ that we have fitted to straight lines, according to eq. ( [49| , i.e. 
without taking into account the finite time extent of the lattice because it affects the correlation functions 
only at large times where no signal can anyway be extracted within our errors. 

In Figure [s] we show the chiral extrapolation of AMjy / Arriud performed by using the following fitting 
function 124| 



AMn 
Arriud 



(rriud) 



AMn 



Arriud 



QCD 



QCD^ 



(50) 



By using the results of the fit and the value of Am^?^ given in eqs. (43), we get 



\QCD 



AM 



N 



Arriud 



QCD 



= 2.8(6) (3) MeV x 



6.05 X 103 MeV^ ' 



(51) 
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where the first error takes into account the lattice uncertainties while the second comes from the uncertainty 
on e^. This number is our best estimate at present but it has been obtained at fixed lattice spacing and 
with limited statistics. We plan to refine the calculation in a separate publication. 



7 Semileptonic decays 

In order to illustrate our calculation of the QCD isospin breaking corrections to the form factors parametriz- 
ing semileptonic decays, we start by the following "double ratio" of three point correlation functions 25 , 



K+Ko it) = ' ' (52) 



where the correlators at numerator are given by 



{ dj^sipi, T/2) s-i^c{pf - Pi, t) C75d(0) ) 
{dj^c(pf,T/2) c^^s{pi-pf,t) 87^^(0)) 



while those at denominator are 



{d^'c{puT/2) C7°c(0,t) C7^d(0)) = {D+\Vi\D+) e'^^-^/^ 



By using the previous expressions and the conservation of fiavour diagonal vector currents, Zy{n\V'^\n) = 
2En^ we obtain the well known and useful result 

K^Ko{t) = ^E^^E^^ \{K^K'^D+)\' + ---, (55) 

i.e. the fact that, by neglecting sub leading exponentials, i?^+^o(^) is a constant with respect to t from 
which it is possible to extract the matrix elements with high statistical accuracy, thanks to the statistical 
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correlation between the different correlation functions. The form factors can be extracted by using the 
standard expressions of the matrix elements computed at different initial and final meson momenta (with 
Pi and/or pf ^ 0) 



(56) 



where as usual ^ (q^) can be expressed in terms of the scalar form factor 



Ml, 



(57) 



Matrix elements with non vanishing momentum transfer between initial and final states have been com- 



puted by using flavour twisted boundary conditions for the valence quarks 26 27 



A very important observation concerning the calculation of isospin breaking corrections is that i?^+^o (^) 
in eq. (55) has the same form, when expressed in terms of diagrams, not only in the perturbed theory (the 
one discussed up to now), but also in the unperturbed isospin symmetric theory. This happens, as in the 
case of two point correlation functions, because C does not generate decays of the D oi K mesons. Thanks 
to this observation we simply have 





25[{K\V,^^\D)]-5Ed-5Ek 



(58) 



From the expressions of the form factors in terms of the matrix elements {K\V^q\D) and by using 1//^ 
and the values of SEjj^k obtained as explained in the previous sections it is possible to extract Sf^^{q^). 



A more detailed derivation of eqs. (58) is provided in Appendix p 



The calculation of 8f^^{q^) proceeds along similar lines but there are some important differences that 
require a separate and detailed discussion. The starting point are the following diagrams for R^j^{t) and 
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its variation, see eqs. (25) and (26), 



1 



4EkE^ 






25[{n\V,^jK)]-5EK 



(59) 



Note the differences of eqs. (59) with respect to eqs. (58), i.e. the presence of disconnected contributions 



and the absence of the correction to the pion correlation function (ah black quark lines) in the denominator 
of Rxni'^)' "^^^ latter is a consequence of the vanishing of the QIB corrections at first order in Arriud in 
the pion case. 

In this work we have not calculated disconnected diagrams and we cannot show results for 8f^^{q^). 
These will be given in a separate publication but, for the time being and in order to show that our 
method works also in the case of three point functions and form factors, we have calculated the difference 
of ^ (q^) with respect the isospin symmetric value f^^{q^)^ i.e 8ff^^{q^). This is a quantity that 
cannot be measured directly because the missing contribution, 5i)f^^{q^)^ is neither equal nor related in 
a simple way to 5ff^^{q^). The two different contributions are in fact associated to two independent 



isospin channels and, according to ref. 28 , the tt^-t] mixing is expected to enhance considerably Si)f^^{q^) 



with respect to Sff^^{q^). One may be tempted, to a first approximation, to neglect the disconnected 
diagrams but this cannot be done because they are needed in order to cancel non physical terms {tAEp) 
contributing to the slope of the connected diagrams. 



= const. - tAEK + tAEp 




= const. - tAEK + • • • . (60) 



The AEp contribution to the slope corresponds to the QIB correction to the energy of a meson, a copy 
of the physical pions, having as valence quarks a physical u (or d) and an additional light quark, also of 
mass niud^ but not contained into the isospin doublet. This term is of the same size of AEk-, as we have 



explicitly checked numerically by using the slopes extracted from two point functions (eq. (30)), whereas it 
cannot be present in physical kaon-to-pion three point correlation functions because of isospin symmetry. 
Indeed the correction to the energy of the physical pions vanishes at first order in Arriud - 
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Figure 6: Left panel: we show our results for y R%j^{t) for several values of the momentum transfer. Right panel: we 

show our results for -^/^^OrCO several values of the momentum transfer. The data are obtained at fixed lattice spacing 
a = 0.085 fm and at fixed am^. = 0.0064 (see Appendix [aI . 
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Figure 7: Left panel: we show our results for dR^j^{t)/aA'm^^ for several values of the momentum transfer. Right panel: 
we show our results for S f R^^{t) / aAm^^ for several values of the momentum transfer. The data are obtained at fixed lattice 
spacing a = 0.085 fm and at fixed am^^ = 0.0064 (see Appendix IaI. 



In order to calculate Sff^^{q^) we need an expression for SfRj^^{t). From eq. (25) we get 




2^/ [{7r\V,^um-SEK 



(61) 



We now come to the numerical results. In Figure [g] we show the ratios R^j^{t), left panel, and 
(t), right panel, from which we extract the form factors. In Figure m we show dR^^j^{t) /aAm^^, left 



panel, and SfR^^{t)/aAm^^. As expected according to eqs. (58) andTGlk both SR%j^{t)/aAm^^ and 



5fR^^{t)/aAm^^ are constant with respect to t in the middle of the lattice, within the statistical errors 
that, in the case of SR^j^{t) /aAm^^^ are rather large. Finally, in Figurejsjwe show our results for f^^{q^)^ 
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Figure 8: Left panel: we show our results for f^'^{q^). Right panel: we show our results for Sff^'^{q^). The data are 
obtained at fixed lattice spacing a = 0.085 fm and at fixed am^^ = 0.0064 (see Appendix [a| . 



left panel, and for 8ff^^{q^). These results, obtained only at fixed lattice spacing a = 0.085 fm and fixed 
light quark mass am^^ = 0.0064, show that our method works also in the case of complicated observables, 
extracted from ratios of integrated three point correlation functions with non vanishing spatial momenta. 

In a separate paper we shall calculate disconnected diagrams and the missing contribution Si)f^^{q^) 
and refine our findings by improving the statistics, by repeating the calculation at several values of rriud 
and by performing chiral and continuum extrapolations. For the time being, we get 



fK . (0)-/K.(0) 

/f-(0) 



QCD 



0.85(18)(1) X 10" 



iQCD 



6.05 X 103 MeV^ 



(62) 



As anticipated a few paragraph above, the fact that this result is two order of magnitude smaller than 
the one quoted in eq. (27) is not surprising in view of the chiral perturbation theory analysis performed 
in ref. ^28j. The enhancement of the missing contribution, ^6/:f"^(<?^), is indeed traced back to the tt^-t] 
mixing. 



8 Conclusions and Outlooks 



In this paper we have proposed a new method to compute with high precision leading QCD isospin breaking 
effects in relevant physical quantities at the lowest non trivial order in the up-down mass difference. 
The method can be easily extended with minor modifications to higher orders. We have computed the 
corrections to meson and nucleon masses, meson decay constants and weak form factors, showing that, 
in spite of the limited statistics, our approach is already competitive, or even better, than other non 
perturbative calculations based on effective chiral lagrangians. 

To obtain the complete physical results, our method has to be combined with calculations of the elec- 
tromagnetic corrections which will be the subject of a future investigation. In this paper, for a comparison 
with calculations in different theoretical frameworks, we have taken the electromagnetic corrections to the 
meson masses evaluated in ref. 
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As the method looks very promising, we are planning to extend this pioneering work to other physical 
observables. 
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Appendix A 



am, 



ud 



am. 



L/a 



conf 



a (fm) 



Zp{MS,2GeV) 



3.80 



0.0080 
0.0110 



0.0194 



24 
24 



150 
150 



0.0977(31) 



0.411(12) 



3.90 



0.0030 
0.0040 
0.0040 
0.0064 
0.0085 
0.0100 



0.0177 



32 
32 
24 
24 
24 
24 



150 
150 
150 
150 
150 
150 



0.0847(23) 



0.437(07) 



4.05 



0.0030 
0.0060 
0.0080 



0.0154 



32 
32 
32 



150 
150 
150 



0.0671(16) 



0.477(06) 



4.20 



0.0020 
0.0065 



0.0129 



48 
32 



100 
150 



0.0536(12) 



0.501(20) 



Table A: Guage ensambels used in this work. The gauge configurations have been generated with Nf = 2 dynamical 
fiavours of maximally twisted quarks of mass am^^. The strange quark mass am J it has been used for valence propagators. 
At (3 = 3.90, where we calculate quantities related to D mesons and correlation functions at non-vanishing spatial momenta, 
we have set am^ = 0.2123 and, by using fiavour twisted boundary conditions, pL/27T = {0.00, ±0.15, ±0.35}. 

In this work we have used the Nf = 2 dynamical gauge ensambles generated and made publicly 
available by the European Twisted Mass Collaboration (see Table [A]) . These gauge configurations have 
been generated by using the so called Twisted Mass lattice discretization of the QCD action [29]. The 
maximally twisted fermion action is given by 



Ctm [U]=q{D[U]+m + ij^T^WlU] ) q , 



(63) 
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where = D[U] is the naive lattice action and W\U] the critical Wilson term (with a = 1), 



V^[C/] q{x) 
V^[C/]t q{x) 

D[U] q{x) 

W[U] q{x) 



U^{x)q{x + Id) - q{x) , 
q{x) -Ul{x- ii)q{x - n) , 



7 



-q{x) , 



E 



v^[;7]-vt[c/] 



q{x) . 



(64) 



The critical mass m^^(^o) been taken from ref. 30 . Concerning our work, the choice of the maximally 
twisted Wilson lattice formulation has advantages and drawbacks. The big advantage is automatic 0(a) 
improvement ^1 . The drawback is the breaking of isospin symmetry at finite lattice spacing even with 
^^lid = 0- Indeed, by letting the physical u and d fields to coincide with the fields ^+ and by taking 
m = ruud and by identifying Co with Cj-m ^ eq. ^ ceases to be valid at fixed cutoff. This happens because 
of the interference between the ts matrix appearing within S and the ts appearing in the twisted critical 
Wilson term. 



[fermionic Wick contractions of O] x tr[/S] ^ 

([fermionic Wick contractions of O] x {Tr [G£^{x,x)] — Tr [G^_(x,x)]} )q 
O(a^) 7^0. (65) 



owing to the invariance of the TM lattice action under 
interchange. Since eq. (65) represent a mere O(a^) cutoff effect we have chosen to 



Notice however that it remains true that {S)o 
parity times ^ £ 

neglect the corresponding contributions to the correlation functions considered in the text. This procedure 
actually corresponds to work within the mixed action approach of ref. 32 and, at the price of introducing 
O(a^) unitarity violations, preserves 0{a) improvement of physical quantities. 



To clarify the point, let's first consider the discretized version of eqs. (22), namely 



Ckoro (t) 



(66) 



In the previous equations we have explicitly shown a label indicating the fiavour of each propagator and the 
corresponding sign of the term ± 75VK[[/] appearing in its kinetic operator. The important point to note is 
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that in the valence we have chosen the signs of the Wilson terms of the up and down quarks independently 
from the choice made in the sea, where the two quarks must necessarily have opposite signs in order to 
deal with a real positive fermionic action, and that this choice is legitimately observable dependent. In 
this particular case the results are correlation functions with much smaller lattice artifacts and statistical 
errors with respect to the other possible choices, e.g. Tr [75GS+ (0, ^)75G^+ (^^ 0)] • 

Analogously we can choose conveniently the sign in front of the Wilson term for the valence quarks 



entering three point correlation functions. As an example we discuss explicitly the case of eq. (25), namely 




"09" 



0{Amud)^ . (67) 



As it can be seen, in order to have mesons interpolated by operators of the form (the ones entailing 

smaller discretization and statistical errors), here we have chosen the sign of the u quark Wilson term 



opposite with respect to that in eqs. (66). The discussion of eq. (26) is considerably more involved because 



of the presence of a neutral pion in the physical correlation functions. A convenient lattice discretization of 



the correlator of eq. ( 26 ) can be also obtained within the mixed action approach but since the calculation 
of hf^^{q^) will be the subject of a future work we don't discuss here this point. 

For the different gauge ensembles used in this work the values of the lattice spacing a (ref. [19'), of 
the strange valence quark mass (ref. |19 ) and of the renormalization constant Zp (ref. 33 ) are given in 
Table [AJ The values of Zp are relevant because in the maximally TM formulation one has [29^^32^ 



Am, 



A qrsq 



m 



1 



(68) 



Appendix B 



A more detailed derivation of eqs. (58), i.e. of 





/\ 




2S[{K\V^^JD)]-5Ed-SEk 



(69) 



can be obtained by applying perturbation theory with respect to Arriud (see refs. 34 36 for related works). 
Let us analyze in detail the case of 



(70) 



25 



where = {^\d'^^ s(fS)\K^) /2Eko and pd+ = (D+|c7^(i(0)|0)/2E^+ and where we have neglected sub 
leading exponentials. The perturbation V, 



V = ^ £(0, x) = Y^ [uu - dd] (0, x) , (71) 



is flavour diagonal and does not open any decay channel for the K and D mesons. Furthermore charged me- 
son states do not mix with the corresponding neutral states. By considering lattice states \nL) normalized 
to one, (tilI^l) = 1, we can use the well known formulae 

0\ _ IZV'X^IAIV'.X _ IIV-.X^A^ ^ \nL){nL\V\KL) 



\Kl) = \Kl) + \^Kl) = |/Ci)+Amw^ 



Eko = Ek + AEk = Ek + Amud{KL\V\KL) , 

\nL){nL\V\DL) 



\D+) = \Dl) + \ADl) = \DL)+AmudYl p_p 
Ed+ = Ed + AEd = Ed + Amud{DL\V\DL) . (72) 

connecting, at first order, the states and the eigenvalues of the unperturbed isospin symmetric theory 
with the corresponding quantities of the perturbed theory. First order corrections for the relativisticaUy 
covariant states, (n|n) = 2£'„, are then readily obtained by changing the normalization. 



= \K) + \AK) = ./2Eko\KI) 

= y^(^l+^y\KL} + \AKL)) 



In the case of \D~^), by proceeding along the same lines, we get 
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The expansion of the matrix element {K^\V^g\D~^) appearing into eq. (70) is thus given by 



(if°IK^|Z)+) 



{K\v,m + + {K\Vc''.\m 



^ {K\V\n) {n\V,m , {K\VrM) {n\V\D) 

2En{En — Ed) 



"^EniEn - Ek) 



E 



(75) 



In order to obtain the exphcit expression of the QIB correction to the correlator of eq. ( 70 ) we should also 



expand the exponentials and the matrix elements of the interpolating operators. The explicit expression 
of SpD^K is not needed, because these terms cancel in the final expression of SR^j^, while the expansion 
of the exponential factors is easily obtained, 



^-tEn^-iT/2-t)EK [1 _ ^^^^ _ - t)AEK] 

{D\V\D) 



^-tED^-{T/2-t)EK 



1 — t Am,, 



2Ek 



(r/2 - t)/\m^ 



{K\V\K) 



2Ek 



(76) 



Eqs. ( [75| and (76), combined with the diagrammatic expansion of our correlation function derived in 
section |3l i.e. 




(77) 



give us an explicit expression for the first term entering the SR^^ formula, namely 




= {1 + 5p),} {1 + 5pd} {1 + 5[{K\V,^^D)]} {1 - t^Eo - (T/2 - t)AEK} - 1 • 



(78) 



By repeating the same arguments for the other correlation functions of eqs. (53) and (54), it is straightfor 



ward to obtain explicit expressions for the remaining three terms appearing into the expression of 5R^ 



DK^ 
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i.e. 



{1 + 5pk} {l + 5pl} {1 + 5[{K\V,^,\D)]} {1 - tAE^ - (T/2 - t)A^c} - 1 ^ 




= {1 + 5pd} {1 + 5p*d} {1 + 5Ed} {1 - T^Ed/2} - 1 , 




{1 + 5pk} {l + 5p)i} {1 + 5Ek} {1 - - 1 . 



(79) 



relation eqs. (78) and eqs. (79). 



The proof of the SR^j^ formula, second of eqs. (58) or eq. (69), is finally obtained by substituting in that 



In deriving eqs. ( |78| ) and eqs. ( |79[ ) we have not shown terms proportional to the correction of the 
"vacuum energy" because in our case (0|V|0) =0. If different from zero, such kind of contributions would 
appear at an intermediate stage of the calculation whereas they would cancel in the final expression as 
happens to the terms proportional to SpD,K- 
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